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Abstract. The contact process is studied on long-range connected networks in the 
presence of random transition rates by means of a strong disorder renormalization 
group method and Monte Carlo simulations. We focus on the case where the connection 
probability decays with the distance x as p{x) ~ Px^"^ in one dimension. Here, 
the graph dimension of the network can be continuously tuned with /3. The critical 
behavior of the model is found to be described by an infinite randomness fixed point 
which manifests itself in logarithmic dynamical scaling. Estimates of the complete 
set of the critical exponents, which are found to vary with the graph dimension, are 
provided by different methods. According to the results, the additional disorder of 
transition rates does not alter the infinite randomness critical behavior induced by 
the disordered topology of the underlying random network. This finding opens up the 
possibility of the application of an alternative tool, the strong disorder renormalization 
group method to dynamical processes on topologically disordered structures. As the 
random transverse-field Ising model falls into the same universality class as the random 
contact process, the results can be directly transferred to that model defined on the 
same networks. 
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1. Introduction 

Dynamical processes occurring on random networks have been in the center of 
interest recently [H |2]. The motivation of this field is many-sided, ranging from the 
spreading of epidemics or rumor in social networks via spreading of computer viruses to 
transportation networks [3] . There are at least two striking properties of the underlying 
networks in these problems which make the dynamics on the top of them much different 
from those on regular lattices. First, the nodes are "close" to each other, meaning that 
the number of nodes n{i) that are reachable by traversing at most £ links from the origin 
is increasing rapidly with i. In the case of an algebraic relationship 

n{i) ~ (1) 

a generalized graph dimension can be associated with the network under consideration. 
In many cases in the above systems the graph dimension of the underlying network is 
formally infinite, like in the case of small world networks [H [5] where. In n{i) ~ i, or 
in scale-free networks with a broad distribution of degrees, where ln[lnn(^)] ~ ^ [3]. 
This "small-worldness" makes mean-field approximations successful in the description 
of phase transitions of spreading processes in many cases [6|. Second, the nodes of the 
underlying network are not equivalent since the degree of nodes, as well as the local 
neighborhoods are different and this circumstance may induce disorder effects on the 
dynamical process occurring on it. This issue has been recently studied [HI Ellin] in the 
case of a paradigmatic model of an epidemic, the contact process In this model, 
there are binary variables attached to each node, which can be either active or inactive 
(infected or healthy in the parlance of epidemiology) . Infected nodes then stochastically 
infect neighboring healthy nodes or recover. Varying the relative rates of these two 
competing processes the system can be driven from an active phase with a finite fraction 
of active nodes to an absorbing one where all sites are inactive. The phase transition on a 
translationally invariant lattice is known to fall into the directed percolation universality 
class [m [m [IE] . This model has been studied on long-range connected (LRC) networks 
which consist of a (i- dimensional hypercubic lattice and additional long-range edges 
existing with probability that decays algebraically with their Euclidean length x as 
p[x) ~ /3a;~*. This type of broad distribution of lengths of links has been observed 
in a mobile phone communication network, where s = 2 [H] and in the global airline 
network, where s = 3 [12]. The extremal cases s = and s = oo correspond to small- 
world networks and short-range networks, respectively. In Refs. [H O |10], where the 
case d = 1 was considered, the appearance of disorder effects was conjectured to be 
related with the finiteness of the graph dimension of the underlying network. This is 
realized in the d = 1 LRC network if the decay exponent s is large enough, s > 2 
namely, which includes the "critical" point s = 2 where the graph dimension varies 
continuously with the prefactor /3 of the asymptotical probability of links [TSj HH |15] . 



I Note that, besides regular lattices, these networks provide an alternative for approaching the limiting 
case dg = oo through tuning the prefactor f5. Nevertheless, the two limiting cases may not be equivalent. 
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In this case, an anomalous slow, algebraic decay of the density of active nodes was 
observed in Monte Carlo simulations of the contact process in an extended phase on the 
subcritical side of the transition point, at least for small enough graph dimension. This 
region is an analogon of Griffiths-McCoy phase of disordered ferromagnets where the 
system, although it is globally paramagnetic, contains locally ferromagnetic domains of 
arbitrary size [19]. In the analogous phase of the contact process on LRC network, the 
majority of the system is locally sub-critical but the randomness of the structure leads 
to the formation of rare regions (sub-graphs) which contain an over-average number 
of internal links and, as a consequence, may be locally super-critical. The activity in 
these rare regions get extinct very slowly, so they give a large contribution to the average 
density and result in anomalous decay. The existence of this phase in the contact process 
with site-dependent random rates on regular lattices has been known for a long time 
[20] . In the following, we will term this type of inhomogeneity as parameter disorder in 
order to distinguish it from topological disorder which refers to the irregularities of the 
underlying network. 

Besides the off-critical behavior, parameter disorder has a striking effect on critical 
scaling, as well [21], where, in the dynamical relations, the time is formally replaced 
by its logarithm. Indeed, an asymptotically exact strong disorder renormalization group 
(SDRG) treatment [23] of the one- dimensional model showed that, at least for sufficiently 
strong disorder, the critical behavior is controlled by an infinite randomness fixed-point 
(IRFP) where the temporal scaling is logarithmic and yielded the complete set of critical 
exponents [22]. Later, this type of critical scaling has been indicated by the SDRG 
method in higher dimensions d = 2,3,4, as well as on the Erdos-Renyi graph [23] 
(where formally dg = oo) with parameter disorder [25l |26l [27]. Based on these results, 
the existence of infinite-randomness critical behavior was conjectured on arbitrarily 
high- dimensional hypercubic lattices for sufficiently strong parameter disorder. 

On LRC networks, where solely topological disorder is present, numerical 
simulations have indicated a logarithmic critical scaling for small enough graph 
dimension, where the critical exponents vary with dg but a precise estimation of them 
is still lacking. Here, the critical point has been found to be flanked with a Griffiths 
phase, the width of which is shrinking with increasing rf^ [Hj |9l [10]. For larger dg 
(from dg ^ 2 on) a Griffiths phase could not be observed and the judgement of the 
critical behavior became less certain. Here the data were compatible with a conventional 
algebraic scaling. 

In this paper, we wish to contribute to the above issue by revisiting some still unclear 
parts, as well as to extend it to another direction. Namely, we shall focus on the critical 
behavior of the contact process on d = 1 LRC networks but, as a new feature, mainly 
in the presence of parameter disorder. Keeping in mind the universality of the IRFP 
in case of parameter disorder, meaning that the critical exponents are independent of 
the distribution random parameters [23], we pose the question whether, on a structure 
with topological disorder, the critical exponents of the transition are modified by an 
additional parameter disorder or not. A methodological yield of this new component 
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is that the model becomes ehgible for a SDRG analysis that always necessitates some 
initial parameter disorder. Furthermore, we expect the additional parameter disorder 
enhancing the effective strength of disorder (i.e. reducing the transient regimes which 
are usually rather long in the IRFP) and making possible a more accurate numerical 
investigation. We shall thus perform a numerical SDRG analysis of the above model 
and compare the predictions with results of Monte Carlo simulations. 

It is worth mentioning that, provided it is controlled by an IRFP, the transition of 
the disordered contact process falls into the same universality class as that of the random 
transverse-field Ising model on any structure, as the strong disorder renormalization 
rules of the two models are identical (apart from factors that are irrelevant at criticality) 
[30] . Besides regular lattices, the (zero-temperature) quantum phase transition in this 
model has been studied also on a LRC network when the underlying network is driven 
through a (long-range) percolation transition [31]. Here, the critical exponents of the 
model can be expressed in terms of those of the percolation transition. However, the 
nature of the quantum phase transition controlled by the strength of the transversal 
field (on the percolating LRC network), which is a much harder problem, and which is 
the subject of this work, has not been revealed yet. Due to universality, the results of 
our investigations can immediately be transferred to the quantum critical behavior of 
the random tranverse-field Ising model. 

The rest of the paper is organized as follows. The precise definition of the model will 
be given in section |2l In section |3l the scaling theory of the IRFP will be recapitulated. 
Section m is devoted to the SDRG analysis of the critical behavior, while results of Monte 
Carlo simulations are presented in section [51 Finally, the results are discussed in section 
|6] and conclusions are drawn in section [7l 

2. The model 

Networks with an algebraically decaying probability of long edges have been studied in 
the past from different aspects. In addition to the examples mentioned so far, they also 
arise as models of linear polymers with crosslinks between remote monomers [32], in the 
context of decentralized search algorithms [33], or, indirectly, in susceptible- infected- 
recovered models with long-range infection [38]. Besides contact process [Sj El [TO] , 
long-range percolation [35], random walks [361 [IS], susceptible-infected-recovered model 
[37] and spin glass models [31] have been studied on them. The geometry of these 
networks itself (the dependence of the diameter on the number of nodes) has attracted 
much interest, as well [131 [391 [iQl [H [HI [121 [M] • 

Concerning the precise definition of the LRC networks, we shall adopt that of Ref. 
[13] . Let us have a set of N nodes, which are labeled by integers 1,2, ... ,N and define 
the distance between node i and j as Xij = min(|i — j\, N — \i — This simply means 
that the nodes are arranged on a ring with unit spacing between them. Then all pairs 
of nodes with a distance Xij = 1 (i.e. neighboring nodes on the ring) are connected with 
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a link and all pairs with Xij > 1 are connected independently with the probability 

p{x) = l-exp{-l3x-'), (2) 

where (3 and s are positive constants. For large x, this probability has the asymptotic 
form 

p{x) ~ (3) 
A finite realization of this network is illustrated in Fig. [H 



Figure 1. Left: A realization of the LRC networks witli N = 256 and 13=1. Riglit: 
The same network represented with roughly equal edge lengths visualized by the Pajek 
program [43j . 

As for the geometry of these networks, for our purposes it is sufficient to know that 
if s = 2, the graph dimension is finite and continously increasing with /3 from one to 
infinity [131 E], for numerical estimates, see Refs. [Ul [3H]- If s > 2, the network is 
short-ranged (i.e. dg = 1), whereas if s < 2, formally dg = oo that hides in the range 
1 < s < 2 an increase of the diameter with a power of IniV In the present work, 
we shall restrict ourselves to the critical value s = 2. 

To the nodes of these networks binary variables coding active {ui = 1) and 
inactive (rij = 0) states are attached, and on this state space a continuous-time stochastic 
process is considered with the following transitions. If node i is active it can become 
inactive with a rate /Xj, as well as it tries to activate each neighboring node j with a rate 
Xij. The rates /ij and Xij are i.i.d. quenched random variables drawn from a distribution 
that will be specified later. 

Note that the homogeneous version of this model, where /ij = /i = const, 
Xij = const = A is known as the susceptible-infected-susceptible (SIS) model rather 
than the contact process. In the former model the infection rates on links whereas in 
the latter the total rate of infection from nodes are constant. In other words, in the 
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homogeneous contact process, the rate of infection from node i through a given hnk is 
X/di, where di is the coordination number (or degree) of node i. Therefore, in the case 
of the contact process on a non-regular lattice the infection rates on a given link can 
be different in the two directions. Although this situation is tractable by the SDRG 
method [34], an efficient algorithm is at our disposal that works for the symmetric case, 
thus the SDRG analysis will be carried out for the SIS model. In computer simulations, 
however, it is more natural to implement the contact process rather than the SIS. (For 
the discrete time- implementation see section El ) Nevertheless, in the case of networks 
with a narrow distribution of degrees, which is the case for the LRC networks with 
an asymptotically A^-independent distribution, the above distinction between the two 
models is irrelevant as far as the critical exponents are concerned. 



3. Scaling theory in an IRFP 

Here, following mainly Ref. [22j, we shall briefly recapitulate the essentials of the scaling 
at an infinite randomness fixed-point, in a form adapted to finite-dimensional networks. 

As it is natural in strongly disordered systems |23] , we define the control parameter 
in the form Aq = ln(A//i), where the over-bar denotes an average over the distribution 
of rates, whereas A = Aq — Aq will denote the deviation from the critical value Aq. As 
we will formulate the scaling theory for random networks, the usual linear size appearing 
in scaling relations will be replaced by the mean diameter D of networks with A^ nodes. 

An appropriate order parameter of the transition is the probability p that a 
randomly chosen node in the stationary state of the infinite system (A^ — )■ oo) is active. 
In the active phase and close to the critical point it vanishes as 

p(A) ~ A^ (4) 

defining the order parameter exponent /3. Close to the transition, the spatial correlation 
length (measured in terms of the Euclidean distance x), ^± diverges as 

a-iAi-'^-, (5) 

with the correlation length exponent u^. The dynamics in an IRFP are strongly 
anisotropic; the time and length scales are, namely, related to each other as 

In^ll-et, (6) 

where i/j is the tunneling exponent and plays the role of a dynamical exponent. Note 
that the conventional dynamical exponent that is defined by an algebraic relationship 
rather than Eq. [6] is formally infinite here. Based on the above, the order parameter for 
finite size D and time t is expected to have the scaling form when the length (measured 
in shortest-path distance £) is rescaled by a factor b as: 

p(D,t, A) = r"p(D/6,lnt/6^, Afoi/"^), (7) 

where x = f3/h'±. 

The observables the time-dependences of which are usually measured in numerical 
simulations starting from a single active node are the survival probability, the average 
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number of active nodes and the spread [IS]. The first one is the probabihty that there 
will be at least one active node at time t. This probability is then averaged over the 
position of the initial active seed. For the contact process, it scales in the same way as 
the order parameter [T71 [18] so we have 

P{D,t,A) = b-''P{D/b,lnt/b^,Ab^/''^). (8) 

Writing a scaling relation for the spatio-temporal correlation function C[no{t = 0), r2j(t)] 
analogous to Eq. ([7]) then summing over the position i yields for the scaling of the 
average number of active nodes: 

Na{D,t,A) = 6'^«-2"iV(D/6,lnt/6'^, Afei/"^). (9) 

The spread is usually defined as the root-mean-square of the Euclidean distance of active 
nodes from the origin. This quantity would, however, diverge in our model for s < 2 
owing to the long links, therefore we define it in terms of the shortest-path distance i 
as 



R{t) = MJ2^.{t)m/J2''^{t)), (10) 

where (■) denotes the expected value conditioned on the survival up to time t for a 
given random environment (i.e. a given random network and given starting position) 
whereas the overbar denotes an average over the latter. Using again the scaling form of 
the spatio-temporal correlation function one obtains that the spread obeys the scaling 
relation 

R{D,t,A) = bR{D/b,\nt/b^,Ab^/''^). (11) 

In the critical point (A = 0) of the infinite system {D = oo) it follows from the 
above relations that the observables depend on time asymptotically as 

P{t) ~ IHt/to)]-', (12) 

Kit) ~ [ln(^/^o)]^ (13) 

Rit) ~ [ln(t/to)]^/^ (14) 

where to is a non-universal microscopic time-scale and the exponents are given in terms 
of the earlier ones as 

5 = x/i), T]={dg- 2x)/i) . (15) 

A further quantity that is directly accessible by the SDRG method is the lowest 
gap e of the rate matrix (or infinitesimal generator) of the process in finite systems (see 
the next section). This quantity varies from sample-to-sample and can be interpreted 
as the inverse of the mean time r needed to reach the absorbing state starting from the 
fully active one in a given finite sample. The distribution of its logarithm in samples 
with fixed iV, in the critical point A = obeys the finite-size-scaling relation 

U\Yit,D) = D-^UD-nnt). (16) 
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Similarly, the distribution of mean lifetime r has the scaling form 

M\nT,D)=D-^UiD-^\nT). (17) 

Note that, in the above scaling relations, we could have equally well used the 
"volume" of networks instead of their diameter D. But in that it is easy 

to see from the relation D(N) ~ N^/'^-a ^ the scaling exponents involving the finite size 
would differ from the present ones by a factor of l/dg and the graph dimension in Eqs. 
flgfTSj) should be replaced by 1. This set of exponents defined in terms of the volume 
will be distinguished by a prime from the above ones: 

'ilj' = ip/dg, x' = x/dg, u' = udg. (is) 

The other exponents are not effected by the choice of the measure of the size. 



4. SDRG treatment 



By the SDRG procedure, the quickly relaxing degrees of freedom - that are related to 
high-lying levels of the rate matrix that governs the time-evolution - are sequentially 
eliminated, while lower-lying levels which are responsible for the long-range dynamics 
are kept |22j . 

To be specific, the procedure consists of two kinds of reduction steps. First, if the 
infection rate Xij on a link is much greater than the recovery rates on nodes i and j then 
the two nodes are merged and treated as a single giant node with the effective recovery 
rate 

ft ~ 2ij,iij,j/Xij (19) 

and size m = nii + nij. Second, if the recovery rate fii on node i is much greater than 
the infection rates on the links emanating from it then node i is eliminated while any 
pairs of nodes neighboring to it are connected by new links with effective infection rates 

Ajfc ^ XjiXik/iii (20) 

on them. 

In the original formulation of the SDRG scheme the largest transition rate is selected 
then either of the above reduction steps are applied and this loop is iterated until the 
system is sufficiently small so that the spectrum of the rate matrix can be directly 
calculated. Although the reduction steps are approximative, they become more and 
more accurate as the renormalization proceeds since the distributions of logarithmic 
rates broaden without limits, and ultimately, the method becomes asymptotically exact 
in the IRFP |30l|23]. 

This renormalization scheme is formally equivalent to that of the random 
transverse-field Ising model defined by the Hamiltonian 

= - 5^ J.X^J - 5^ /..af , (21) 
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where a^'^ are Pauli operators on site i, Jij and hi are random couplings and external 
fields, respectively, and the first sum goes over neighboring nodes of the network [30] . 
The asymptotic renormalization rules of the model take the forms given in Eqs. ( 1T9|20|1 
with the correspondences /i, -H- hi and Xij O Jij, apart from the absence of the factor 
of 2 in Eq. (fT9l) . which, however, does not influence the critical exponents at an IRFP. 

The performance of the SDRG method in the above form may be rather ineffective 
in other than one dimension since the connectedness of the underlying network will 
rapidly increase by the renormalization. To avoid this we will use a more efficient 
algorithm developed by one of us [261 12H], which is based on the so called 'maximum 
rule'. According to this, if multiple links between nodes would appear during the 
renormalization, only the one with the maximal rate is kept. Application of this rule is 
thought not to influence the critical exponents at an IRFP and results in simplifications 
of the SDRG procedure. The improved algorithm works by merely deleting links (and 
changing the rates on the remaining ones), without generating new ones at all. The 
results of the algorithm are identical to that of any traditional implementation of the 
SDRG method (having also the widely applied maximum rule) for any finite graphs, 
with sites and E edges. However, we gain considerable time in performance: while 
the traditional method needs 0{N^) operations, the improved algorithm requires only 
0{N log N + E), which is in practice much faster, than the (9(A^^) operations needed 
to generate the LRC networks. The essence of the algorithm is that if two (or more) 
sources are able to infect each other mutually before being healthy again, than these 
form a new effective infection source, for which it takes a longer time to become healthy 
again. So, in this sense, the process can be regarded as a special kind of percolation of 
the infection sources with a positive feedback [26l |28] . 

In the numerical SDRG analysis, both variables Ajj and fii where taken from 
uniform distributions with probability densities /a(A) = 6(A)6(1 — A) and /^(yu) = 
— 0(//)0(/Xm — f'), respectively, where 6(x) is the Heaviside step-function and we used 
the logarithmic variable, Aq = ln(/im) as a control parameter. 

It is generally a challenging task to precisely locate the critical point, and the 
accuracy of the obtained critical exponents depends crucially on it. In order to reach a 
sufficient precision, we have first determined the location of the pseudo-critical point 
for each random sample, where the correlation length reaches the size of the system. 
This can be conveniently obtained by the doubling method [291 [28]. Here, the originally 
periodic system is made "open" by cutting all links going over a given position x, and 
two identical copies of this open system are glued together at the boundaries. Then, A^ 
is given by the threshold value, above which the last remaining giant nodes (or clusters) 
of the replicas fuse together during the SDRG process. A^ varies from sample to sample, 
but from the size-dependence of its distribution both the location Aq of the 'true' critical 
point and the correlation length exponent i>' can be obtained according to the scaling 
form: 

/A(A,Ar)=ArV-7^(AAri/'^'), (22) 
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where A = Ac — Aq. In practice, we study the scahng of the width of the distribution 
(given by the standard deviation), which is proportional to N~^^'^^'^ and the mean value, 
|Aq — Ac| ~ N~^^^!''^, and calculate size- dependent effective exponents by two-point and 
three-point fits, which are then extrapolated as — )■ 0, see Fig. [2l 



1.2 




° 10 100 1^ 1000 10000 



Figure 2. The distribution of the pseudo-critical points at /3 = 0.75 is illustrated by 
the mean value and the width of the distribution for various values of the system size 
N. The extrapolated value of the critical point is Aq — 1.1606(15), indicated by the 
solid line. Inset: The obtained finite-size estimates of the exponents 1/i^' lead to the 
asymptotic value 1/v' = 0.308(8). 

Having at hand an accurate estimate for the location of the critical point, the 
remaining two independent critical exponents {x' and can be determined by applying 
the SDRG method at the critical point and analyzing the resulting cluster structure. 
The gap e of the rate matrix (or energy gap in the language of the Ising model) is given 
by the effective recovery rate of the the last decimated giant node (cluster). The mass 
m of the latter gives the magnetic moment of the sample in the corresponding Ising 
model, for details, see [25l l26l |28] . The number of independent random samples used in 
the calculations were at least 40000, while the largest analyzed system size was typically 
N = 2^^ 

According to the scaling theory presented in the previos section, the average mass 
scales with as m(N) ~ N^~^' , since it is related to the order parameter (magnetization 
of the Ising model) p as p = m/N. This is illustrated in Figj3] for /3 = 0.75. For 
sufficiently large system sizes, the points fit well to a straight line. From two-point fits 
we can obtain finite-size estimates of 1 — x', which is then extrapolated as — > 0. 

The mean value and the width of the distribution of In e is shown in FigJH for 
P = 0.75. As a clear indication of infinite disorder scaling, the width of the distribution 
is increasing with A^. According to the scaling form in Eq. [161 both the mean value and 
the width is asymptotically proportional to A^'^'. Similarly to the other exponents, we 
determined finite-size estimates for ip' through two-point fits, as illustrated in the inset 
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N 



Figure 3. The average mass of the last cluster as the function of the system size, N 
at /3 = 0.75. Inset: The obtained finite-size estimates of 1 — lead to the asymptotic 
value 1-x' = 0.575(10). 




100 ^ 1000 10000 



Figure 4. The mean value and width of the distribution of the logarithmic gap plotted 
against the system size N for /S — 0.75. Inset: The obtained finite-size estimates of ip' 
lead to the asymptotic value ip' — 0.32(3). 

of Fig. m The estimates obtained by extrapolations to iV — oo are presented in Table 

m 

5. Monte Carlo simulations 

We have performed Monte Carlo simulations of the contact process on LRC networks, 
implemented in the following way. A finite LRC network with N nodes with i.i.d. 
random variables Wij = Wji on each link have been generated. Choosing an active node 
(z) randomly, it is set to inactive with a probability 1/(1 -|- A). Otherwise, a neighboring 
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ih 


0.1 


1.104(2) 


0.1738(4) 


0.422(7) 


0.214(3) 


0.524(30) 


0.58(3) 


0.2 


1.212(4) 


0.3505(5) 


0.366(4) 


0.231(5) 


0.53(3) 


0.64(4) 


0.3 


1.353(7) 


0.529(2) 


0.323(10) 


0.264(6) 


0.525(20) 


0.71(3) 


0.4 


1.499(7) 


0.700(1) 


0.301(6) 


0.291(10) 


0.485(30) 


0.73(5) 


0.5 


1.656(8) 


0.856(1) 


0.295(8) 


0.329(7) 


0.448(35) 


0.74(6) 


0.75 


2.03(2) 


1.1606(15) 


0.308(8) 


0.426(8) 


0.33(3) 


0.67(6) 


1.0 


2.347(17) 


1.3685(10) 


0.333(10) 


0.469(6) 


0.277(17) 


0.65(4) 


2.0 


3.045(27) 


1.853(3) 


0.335(4) 


0.531(10) 


0.22(4) 


0.67(12) 



Table 1. Critical exponents estimated by the SDRG method for different values of /3. 
The graph dimensions are taken from Ref. [15j; for more precise estimates, see Ref. 

m- 

node (j) is chosen equiprobably and it is activated with a probabihty Wij (provided it 
was previously inactive). The random variables Wij have been drawn from a discrete, 
binary distribution with probability density f{w) = c5{w — wq) + (1 — c)5{w — 1). 
This can be interpreted in a way that a fraction c of the links has a reduced capacity 
of transmitting the disease that is characterized by the parameter wq < 1. The 
parameters of this distribution have been chosen to be c = 0.5 and wq = 0.2 throughout 
the numerical simulations, except of the homogeneous model (i.e. the model without 
parameter disorder) where formally c = 0. One Monte Carlo step of unit time consists 
of Na(t) such updates where Na(t) is the number of active nodes at the beginning of the 
step. We have generated random networks with a fixed /3 in the range [0.2, 2] and with 
= 1 — 5 ■ 10^ nodes (larger ones for larger dg) and - starting with a single active seed 
- we have simulated the process and measured the survival probability, the number of 
active nodes and the spread as a function of time, typically, for MC times up to 2^°. 
This measurement has been repeated for a fixed (3 and A in 10^ independent random 
networks, starting the process from 10^ different nodes per sample, and the measured 
data have been averaged. When presenting results of the Monte Carlo simulations, we 
will simply use A as a control parameter rather than Aq defined in section [31 

The critical point have been estimated in the way proposed in Ref. pEj in order to 
avoid the difficulties about the large time scale to in Eqs. (fT3!) . Plotting In Na(t) against 
InP(t), the slope of the curve must tend to +1 in the inactive phase, to —f]/6 in the 
critical point, and to — oo in the active phase, see Fig. [51 Having identified the critical 
point for different (3, we have determined fj/S hj a linear fit to the critical curve, see 
Fig. [HI Using Eqs. f lT5|TS]) . the exponent x' has then been calculated from 

x' = l/(2 + r//5). (23) 

The estimates obtained this way are given in Table [H Note that the most significant 
source of the error comes from the uncertainty of the critical point. 

The tunneling exponent ijj has been determined from the dependence of the spread 




Figure 5. The logarithm of the average number of active nodes plotted against 
the logarithm of the survival probability for different values of A for (3 = 0.3 (left) 
and (3 = 1 (right). The critical point is located at A* 5.22(1) for f3 = 0.5 and at 
A* = 3.408(3) for /3 = 1. 




Figure 6. The ratio of exponents rj/S measured in Monte Carlo simulations for 
different values of /?. The value at /3 = is the SDRG prediction for the one-dimensional 
random contact process [22] . 



on time by fitting a function R{t) = a[lio.{t/tQ)Y^'^ + 6 to tlie numerical data for times 
\nt > 5. For an illustration, see Fig. [71 The estimates for various /3 can be found in 
Table O For /3 = 1, we have determined the tunneling exponent also by measuring the 
lifetime r starting from a fully active initial state in finite systems. Constructing the 
histogram of In r for different values of A^, the exponent ip can be determined, according 
to Eq. ( 1T71) . by finding an optimal scaling collapse of the data, see Fig. [71 The obtained 
estimate {ip' = 0.23) is compatible with that obtained from the time-dependence of the 
spread. As can be seen in Fig. [HI the dependence of the survival probability and the 
average number of active nodes on time is in agreement with the logarithmic scaling 
laws given in Eqs. ( [T211T31) . 

In addition to this, we have performed simulations and measurements for model 
without parameter disorder (i.e. formally c = 0), as well, in the same way as for the 
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Table 2. Critical exponents estimated by Monte Carlo simulations for different values 
of (3. The graph dimensions are taken from Ref. [T3]. Estimates in the absence of 
parameter disorder (c = 0) are shown in the last three columns. 




[ln(t/to)]''^ ^ ln(t/To)N-^'' 



Figure 7. Left; Time-dependence of the spread measured in numerical simulations 
for P — 1, a.t X — 3.410. The solid line is a fit to the data, yielding the parameters 
Into = 2.4(4) and ip = 0.53(4). Right: Scaling collapse of the histogram of logarithmic 
lifetimes obtained for different sizes N. For each N, the lifetime has been measured in 
10"* — 10^ random samples, once in each sample. Optimal collapse is achieved by the 
parameters Inio = 3, ■(/'' = 0.23. The inset shows the unsealed histograms. 

model with parameter disorder. The obtained estimates are shown in the last three 
columns of Table El 

6. Discussion 

We have studied the critical contact process and the slightly different SIS model in 
the presence of parameter disorder on finite- dimensional random LRC networks by an 
SDRG method and by Monte Carlo simulations. We have seen by both methods that 
the behavior of different observables are compatible with the scaling laws in infinite 
randomness fixed points and have obtained estimates of the critical exponents. As can 
be seen in Fig. |9l they vary smoothly with the graph dimension of networks. We have 
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30 




Figure 8. Time-dependence of tlie survival probability and the average number of 
active nodes for /3 — 1. According to Eqs. (jl2ll3[) . the asymptote of the critical curve 
must be a line. The estimates of the critical exponents 5 = 1.94, rj — 0.54 have been 
calculated from the data in Table [2] using Eq. ([TS]) . 




observed the infinite-randomness scenario all the way to the largest value of dg considered 
in this work and, taking into account the smooth change of critical exponents with dg 
we conjecture this scenario to persist to arbitrarily large, finite graph dimensions. 

In addition to the model with parameter disorder we performed simulations of the 
model with homogeneous parameters. Comparing the estimates of the critical exponents 
obtained in the two cases (see Table [2]) we can see that they agree within the error of 
measurements. This suggests that, at least in the model under study, the presence of 
an additional parameter disorder does not alter the critical exponents of the infinite 
randomness fixed point induced solely by the topological disorder. For supporting this 
finding, a non-rigorous argument based on the SDRG approach of the problem can be 
provided, as follows. 

Let us consider the SIS model and the contact process without parameter disorder. 
All recovery rates are now equal and, in the former, even the infection rates are uniform. 
The SDRG procedure, as formulated for parameter disorder, investigates small blocks 



Infinite randomness critical behavior of the contact process 



16 



(consisting of two nodes in case of the fusion step) and tries to replace them by a 
single effective node. The success of the method lies (among others) in that a local 
control parameter, which is determined by the relative magnitude of local recovery and 
infection rates within the blocks exists. This is much different for topological disorder. 
Here, the locally supercritical regions are expected to be those sub-graphs which have 
an over-average number of internal links. These regions cannot be detected simply 
by investigating and reducing two-node blocks sequentially in the way the procedure 
works for parameter disorder. Instead larger blocks (at least three nodes) should be 
investigated so that the clusters which are super-critical due to the large number of 
internal links can be found Such a hypothetical SDRG method would be certainly 
much more complicated than the naive one but the difficulties are merely of technical 
nature and not principal ones. So we can assume that a renormalization procedure 
exists by which richly linked sub-graphs are selected and replaced by a single node 
with a properly calculated effective recovery rate. Obviously, the application of this 
pre-renormalization leads to that parameters of the model become disordered. In other 
words, we expect the topological disorder of the original model inducing parameter 
disorder in the coarse-grained (renormalized) model. Performing the hypothetical 
SDRG method up to a finite length scale, the resulting coarse-grained model will 
therefore be suitable for an SDRG analysis by the naive method developed for parameter 
disorder. Considering that the pre-renormalization carried out up to a finite-scale does 
not infiuence the asymptotical properties of the geometry of the underlying network 
and keeping in mind the universality in an IRFP, we conclude that the additional 
parameter disorder will not alter the critical behavior with respect to those induced 
by the topological disorder. 

Let us now return to the comparison of the SDRG and Monte Carlo results. 
Concerning the exponent x', the data obtained by the two methods are close to each 
other for small dg but the deviation is significant and increasing for larger dg, see Fig. 
m However, for larger dg, the effective strength of disorder is weaker, which manifests 
itself in that ip'{dg) is decreasing with dg, and this makes the estimates by both methods 
less reliable for larger graph dimensions. Regarding that the agreement is quite good 
for moderate dg, we conjecture that both models are in the same universality class for 
any dg and attribute the discrepancy for larger dg to possible systematic errors related 
to the finite size of networks used in the numerical analysis [Jj]. As can be seen in Fig. 
[9l the exponents x' and u' vary with the graph dimension slowly for large dg and the 
exponents u' fulfill the rigorous bound u' > 2 for all dg |17j. The dependence of these 
exponents on dg is similar to that on the dimension in hypercubic lattices [25]. 

The tunneling exponent ip is found to saturate for large graph dimensions according 
to both methods, see Table [1] and [21 although the values obtained by the SDRG 

§ This problem also arises in the case of the contact process on diluted hypercubic lattices. Dilution 
and parameter disorder have been observed to be in the same universality class 1451 . 
II Note that deviations between the SDRG and Monte Carlo estimates have been established also in 
the case of the contact process on hypercubic lattices with parameter disorder [46] . 
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method are systematically higher than those obtained by Monte Carlo simulations. 
For the exponent ip, one can easily establish an upper bound for arbitrary /3, as follows. 
According to Eq. ([Hj), the exponent ijj' governs the scaling of the typical time t during 
which the activity spreads (in surviving trials) in a finite sample of size N from the 
origin to the node in the Euclidean distance N/2, through Int ~ N^' . Furthermore, 
increasing /3 amounts to that the number of long links increases in the network. Then, 
it is plausible to assume that the time t must not increase with increasing /3 since the 
long edges promote the spreading of activity. Consequently, the exponent ip'{(3) must 
not increase with /3. Since = 0) = 1/2, it follows that 

< 1/2 (24) 

for any /3. 

As can be seen from the data, the decreasing tendency of and the above 

inequality are fulfilled except of the SDRG estimates for small values of /3; nevertheless, 
the upper bound lies within the error of estimates also here. This may be attributed to 
the relatively strong finite-size corrections, depending heavily on the chosen distribution 
of the parameters. In order to have more accurate results, even larger system sizes should 
be needed, or, possibly, other forms of disorder, where the finite-size corrections have a 
different sign. 

It is interesting to note that the possible saturation of ip with increasing dimension 
has been observed also on hypercubic lattices, although to a value (0.46) lower than 
that obtained by either methods of this work on LRC networks |25j . 

The exponents describing the critical dynamics of the contact process, S and fj, can 
be calculated from the data in Table [T] and [H Both vary with the graph dimension, the 
former increasing, the latter decreasing with dg, in agreement with the tendencies on 
hypercubic lattices. 

7. Conclusions and outlook 

We have studied in this work the contact process, as well as the SIS model on LRC 
networks, where the topological disorder leads to infinite-randomness critical behavior. 
We have found by a numerical analysis supported by heuristic argumentations, that an 
additional parameter disorder does not change the critical exponents of the transition. 
In this way, the exponents are universal in the sense that they are determined exclusively 
by the topology of the network. This result opens up an alternative way of investigating 
infinite randomness critical behavior induced by topological disorder in general. Namely, 
after introducing parameter disorder, which is irrelevant in the above sense, one can 
apply the efficient SDRG method to the model. 

Furthermore, the present numerical study has shown that in the case it is 
questionable whether the phase transition in a given model is of IRFP type or 
conventional one, it is better to concentrate on the critical behavior rather than searching 
for a Griffiths phase (an accompanion of IRFP), which may be, in case of a relatively 
weak disorder (small ip'), very hard to detect. 
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We have clearly seen in the present work, that the critical exponents smoothly vary 
with the topological dimension and, for moderate values of dg we provided estimates on 
them. We have obtained from these data indications on the possible limiting behavior 
of these exponents when — )■ oo. Although the exponents for finite dg does not seem 
to be determined exclusively by dg (cf. the estimates on hypercubic lattices [25]), it is 
an intriguing question whether the limiting values of them are universal for dg ^ oo. 
In a wider context, it is also a challenging question, whether there are such topological 
characteristics of the underlying networks which determine the critical exponent of the 
IRFP on the top of them unambiguously or at least approximately. 

As mentioned in the Introduction, the transverse-field Ising model is in the same 
universality class as the contact process, in the presence of parameter disorder. Our 
results indicate that the critical behavior of the above model defined on LRC networks 
where there is exclusively topological disorder would be also described by an IRFP 
characterized by critical exponents obtained in the present work. 
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